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Abstract 

The leading power asymptotic behaviour of the dimensionally regu- 
larized massless on-shell planar triple box diagram in the Regge limit 
t/s — > is analytically evaluated. 
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Systematical analytical evaluation of two-loop Feynman diagrams with four ex- 
ternal lines within dimensional regularization |l| began three years ago. In the pure 
massless case with all end-points on-shell, i.e. pf = 0, % = 1,2,3,4, the problem of 
analytical evaluation of two-loop four-point diagrams in expansion in e = (4 — d)/2, 
where d is the space-time dimension, has been completely solved in |2|, ||, |], |], [6|, |7j . 
The corresponding analytical algorithms have been successfully applied to the evalua- 
tion of two-loop virtual corrections to various scattering processes || in the zero-mass 
approximation. 

In the case of massless two-loop four-point diagrams with one leg off-shell the 



problem of the evaluation has been solved in R 10 , with subsequent applications 



[|TTJ to the process e + e~ — > 3jets. (See |12[ for recent reviews of the present status of 
NNLO calculations. See (13| for a brief review of results on the analytical evaluation of 
various double-box Feynman integrals and the corresponding methods of evaluation.) 
For another three-scale calculational problem, where all four legs are on-shell and 
there is a non-zero internal mass, a first analytical result was obtained in |L4| for the 
scalar master double box. 

The purpose of this paper is to turn attention to three-loop on-shell massless 
four-point diagrams. As a first step, the leading power asymptotic behaviour of the 
dimensionally regularized massless on-shell planar triple box diagram shown in Fig. 1 
in the Regge limit t/s — > will be analytically evaluated. This calculation will 
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Figure 1: Planar triple box diagram. 

demonstrate that a three-loop BFKL analysis [[HJ (at least its virtual part which can 
be reduced to the evaluation of Regge asymptotics) is possible^. 

The calculation will be based on the technique of alpha parameters and Mellin- 
Barnes (MB) representation which was successfully used in [TJ, || |[| |TJ| and reduces, 
due to taking residues and shifting contours, to a decomposition of a given MB integral 
into pieces where a Laurent expansion of the integrand in e becomes possible. At a 
final stage, summation formulae for series of Si(n)S3(n)/n 2 , ip"(n + l)S2(n)/n etc., 
where Sk{n) = Y^=\j~ k i are used. A table of such formulae is presented in Appendix. 

2 In three loops, non-planar diagrams as well as higher terms of expansion of double boxes in e 
are also needed. 
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The general planar triple box Feynman integral without numerator takes the form 
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where s = (pi + P2) 2 and t = (p 2 + P3) 2 are Mandelstam variables, and fc, I and r are 
loop momenta. Usual prescriptions k 2 = k 2 + iO, s = s + z'O, etc. are implied. 

To evaluate the leading power asymptotic behaviour of the master triple box flU), 
i.e. for all ai = 1, in the limit t/s — > one can use the strategy of expansion by 



regions [16|, |17[ 18]. It shows that in the leading power only (lc-lc-lc) and (2c-2c-2c) 



regions contribute, with the leading power behaviour 1/t. (See JIT] and Chapter 8 of 



18 for definitions of these contributions.) The leading power (2c-2c-2c) contribution 
the general planar triple bo> 

T^- 2c - 2c \a u ...,a l0 ;s } t;e) 



for the general planar triple box takes the form 
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where p is such that p 2 = t, 1p\p = 0, 2p 2 p = —t. The leading power (lc-lc-lc) 
contribution is obtained due to the symmetry {1 <-> 3, 4^6, 8^ 10}. 

To resolve the singularity structure of Feynman integrals in e it is very useful to 
apply the MB representation 

1 11 r+ioo Y z 

dz—-r(v + z)T(-z), (3) 



(X + Y) v T{v) 2m J-ioo X"+ 

that makes it possible to replace sums of terms raised to some power by their products 
in some powers, at the cost of introducing extra integrations. In 0, [|, || MB inte- 
grations were introduced directly in alpha/Feynman parametric integrals. It turns 
out more convenient to follow (as in || [L4|]) the strategy of [|19j and introduce, in a 
suitable way, MB integrations, first, after integration over one of the loop momenta, 
r, then after integration over I, and complete this procedure after integration over 
the loop momentum, k. 

After appropriate changes of variables we arrive at the following fourfold MB 
representation of (@): 

(i7T d/2 Y (-l) a 



T(4 - a 589( 10) - 2e)(-s) a 348(_t)«125679(10)-6+3e 



2 



llj=2, 5, 7,8,9,10 1 l a jJ \l1Xl) J-tco j=1 
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where a = 52i=i a i; ^589(10) = a 5 + a8 + a 9 + aio, ct3 4 8 = fl3 + a 4 + a8, etc., and integration 
contours are chosen in the standard way. 

We then turn to the master triple box, i.e. set = 1 and apply, after a preliminary 
analysis, a standard procedure of taking residues and shifting contours which results 
in contributions where one can expand in e in the integrand. In this way, it is not 
immediately clear how to perform integrations in some of 2-dimensional MB integrals 
obtained. One can however proceed numerically at this stage and arrive at a result 
in expansion in e where the simple pole and the finite part are numerically evaluated. 

On the other hand, one can follow a straightforward method based on the MB 
representation and organize the evaluation of the leading power asymptotic behaviour 
in such a way that it will be a part of a future evaluation of the unexpanded triple 
box. Using the same way of introducing MB integrations as outlined above, after 
appropriate changes of variables, we arrive at the following (only!) sevenfold MB 
representation of the (unexpanded) triple box (|l|): 
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In the case of the master triple box, we set a* = 1 for i = 1, 2, . . . , 10 and obtain 

r< >(s,t;e)=r(l,...,l;s,t;e) 
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(in d/2 Y 1 /•+**> * fty r(l + w)T{-w) 
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r(i - z 5 )r(i - Z6)r(i - 2e + z 5 + ^ + *?) 
x r(-i - e - z 5 - z 7 )r(-i - e - z 6 - z r )r(l + z 7 ) 

xT(l + e + w - z 4 - z 5 - z 6 - z 7 )V(—e - w - z 2 + z 5 + z 7 ) 

xT(-e - w - z 3 + z 6 + z 7 )T(l + z 5 + z 6 + z 7 )T(2 + e + z 5 + z 6 + z 7 ) . (6) 

Observe that, because of the presence of the factor r(— 2e) in the denominator, we 
are forced to take some residue in order to arrive at a non-zero result at e = 0, so 
that the integral is effectively sixfold. 

The asymptotic expansion in the Regge limit is determined, in the language of 
the MB representation (Q), by the poles in the variable w that come from a —w 
dependence of the gamma functions. The poles of F(—w) correspond to the hard 
contribution to the asymptotic expansion (within expansion by regions) which turns 
out to be subleading (order t°). The rest of the poles are not explicit and are generated 
due to integration over the other variables, z%. An analysis of the integrand shows that 
the leading Regge asymptotic behaviour, which is effectively described by a gamma 
function of the type T(— 1 — 3e — w), is generated by the gamma functions 

r(— e + z 2 )T(-e - w - z 2 + z 6 + z 7 )T(-1 — e — z% — z 7 ) 

or, symmetrically, by 

r(-e + z 3 )T(-e -w- z 3 + z 5 + z 7 )T(-l - e - z 5 - z 7 ). 

Then the standard procedure of shifting contours and taking residues is applied. It 
results again in a sum of MB integrals where a Laurent expansion in e in the integrand 
becomes possible. 

Following this second variant of evaluation we see that all integration over Zi 
variables are then explicitly performed by means of the first and the second Barnes 
lemmas and their corollaries. The last integral, over w, is taken by closing the inte- 
gration contour in the w complex plane to the right and taking a series of residues at 
w = 0, 1, 2, ... At the final stage, summation formulae for series involving 1/n 3 , Sk(n), 
Sikip) = Yll=ij~ l Sk(j) i etc- are used. Some of them are presented in Appendix. 

The final result for the Regge asymptotics of the planar triple box takes the 
following form: 
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(7) 



where L = lns/t and ((z) is the Riemann zeta function. 

A non-trivial check of this result comes from the first variant of the evaluation of 
the Regge asymptotic presented above where the coefficients at the poles l/e- 7 , j = 
2, . . . , 6 agree exactly with (0) while the coefficient at 1/e and the finite part in e 
agree numerically. 

The value 16/9 for the coefficient at the highest pole is also in agreement with the 
explicit result for a general L-loop ladder planar diagram^ (L-box) by Davydychev 
and Ussyukina [B5 if we believe in a well-known empirical rule which corresponds 
1/e to a logarithm of an expansion parameter. From their result one can find the 
leading power and logarithmic asymptotics: (m 2 ) L A L /(L\s L t) ln 2L (— s). This gives, 
in particular, 4 for L = 1 and L = 2, and 16/9 for L = 3, in agreement with 
existing results at one-, two- and (now) at three- loop level. For higher poles, this 
correspondence hardly exists, even in some generalized form. 



The coefficients at 1/e 6 ,..., 1/e 3 in 



jj) have been also confirmed by a numerical 
based on numerical integration in the space 



check [20| with the help of a method 
of alpha parameters. 

The procedure described above can be applied, in a similar way, to the calculation 
of Regge asymptotics of any massless planar on-shell triple box. 
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00-00313, and the Volkswagen Foundation, contract No. 1/77788. 



A Summation Formulae 

Summation formulae with at least 1/n 2 factor, which explicitly provides convergence, 
for example, 

J2Si(n- l)S 12 (n - l)^ = - 2C(3) 2 , 

~^ n z 45360 

can be obtained^] using FORM (22| procedures called SUMMER and described in ||23|| . 

3 This result stays unique, starting from the two-loop level, for the class of massless four-point 
diagrams with all legs off-shell. For example, if one contracts some line (other than a rung) or puts 
a dot on some line, no analytical results are available for the moment. 

4 I have derived these formulae myself and then was informed about existence of these procedures. 
For a person who has never used FORM, it looks simpler to derive (and check) these formulae than 
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Here is a table of summation formulae with the factor 1/n. They are not explicitly 
present in SUMMER. The convergence is here provided by other factors, ip^(n+l) = 
{-l) k k\ [S k+1 (n) - ((k + 1)] , A; = 1,2,.... We have 
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